o 
o 

(N 
X> 

PL, 

in 



43 

o 



The Role of the Molecular Association in the Change of the Sign of the Soret 

Coefficient in Aqueous Mixtures 

Bernard RousseatQ 

Laboratoire de Chimie-Physique, Bailment 350, Universite Paris-Sud, 91405 Orsay Cedex, France 



Carlos Nieto-Draghi and Josep Bonet Avalos 
Departament d'Enginyeria Qmmica, ETSEQ, Universitat Rovira i 
Avda. dels Paisos Catalans 26, 43007 Tarragona Spam 
(Dated: February 2, 2008) 



Virgili, 



The change of sign with composition in aqueous mixtures of associating fluids is analysed by 
means of Molecular Dynamics simulations. The results obtained are in quantitative agreement with 
the experimental data in water-ethanol and water-methanol solutions, which exhibit the mentioned 
change of sign. A subsequent theoretical analysis is addressed to establish a relationship between 
the dependence of the Soret coefficient with composition with the existence of large inter-species 
interactions. Although the change of sign of the Soret coefficient with composition is not only due 
to so-called chemical contribution analysed here, we discuss the role that these interactions play in 
such a change of sign. 
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Aqueous solutions often exhibit a non-ideal behaviour 
in both thermodynamic and in dynamic properties, due 
to the existence of strong, very directional hydrogen- 
bond interactions 0, Q • One of the most astonishing fea- 
tures observed in many aqueous mixtures is the change of 
the sign of the thermal diffusion factor, or Ludwig-Soret 
coefficient, with composition This phenomenon has 
been observed since long time ago but despite many ef- 
forts there is not even a qualitative microscopic picture 
of the effect yet. The work addressed in this letter is 
twofold: on one hand, we present, to the best of our 
knowledge, the first molecular dynamics simulations able 
to quantitatively reproduce the experimental values of 
the thermal diffusion factor in aqueous solutions of asso- 
ciating fluids, including the change of sign of this coeffi- 
cient with composition. In the light of these results and 
MD simulations of mixtures of Lennard- Jones particules 
of hypothetical compounds, we analyse the role played by 
strong cross-interactions on the composition dependence 
of the thermal diffusion coefficient, and its influence on 
its change of sign. 

When a binary mixture is under a thermal gradient, 
one component will enrich in the cold region while the 
other will migrate towards the hot boundary. At the sta- 
tionnary state, the magnitude of the separation is char- 
acterised by the thermal diffusion factor 



T 



ax 



xa{1 — xa) VT 



(1) 



where xa is the molar fraction of species A. The Soret 
coefficient, St = Q-t/T is also commonly used to quan- 
tify the amplitude of the separation. The sign of ax has a 
very important meaning: it reflects the direction of sepa- 
ration of the components of the mixture when a thermal 
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gradient is applied. The generally adopted convention 
is to take component A as the heaviest one. Thus, ar 
is positive when the heaviest species (A) enrich at the 
cold side, i.e. Vxa and VT have opposite sign. In gas 
mixtures, cut is proportional to the mass difference be- 
tween species, is positive, although this depends on the 
interaction potential Q , and mostly independent of com- 
position. The situation is less simple in dense fluids. In 
a recent paper, Debuschewitz and Kohler [5j have shown 
that the Soret coefficient can be splitted into indepen- 
dent contributions: a mass effect, depending on the mass 
and moment of inertia difference of the species, and a 
so-called chemical effect. The mass effect is found to 
be independent of composition, while the composition 
effect comes solely from the chemical part. In benzene- 
cyclohexane isotopic mixtures, a change of sign was even 
observed at a given composition. In mixtures of associ- 
ating fluids, in some ternary polymer solutions 0,0 and 
electrolyte solutions §, ctT strongly depends on com- 
position. Therefore, the composition dependence of the 
thermal diffusion factor seems to be completely controled 
by the different molecular interactions between the con- 
stituants, which has been referred to as chemical effect. 

To analyse this problem in more depth, we have per- 
formed molecular dynamics simulations on several mix- 
tures where strong interactions exist. We present studies 
of water-ethanol and water-methanol mixtures for which 
experimental data are available. Moreover, we have also 
investigated water-acetone and water-dimethyl-sulfoxide 
mixtures, for which we give our predictions. For all 
studied systems, a strong composition dependence and 
a change of sign is observed. 

Previous applications of MD simulations to the study 
of thermal diffusion have been mostly devoted to the de- 
velopment of methodologies IsL JToL Jill and to the study 
of Lennard- Jones like mixtures [12lll3| . There have been 
only a few attempts to compute the thermal diffusion fac- 
tor in molecular liquids [TH Il5l]. Recently, however, the 
ability of MD simulations in quantitatively reproducing 
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thermal diffusion in pentane-n-decane mixtures has been 
confirmed 
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an important dependence of St with composition is ob- 
served in these mixtures having strong hydrogen-bonds 
interactions. 
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FIG. 1: Thermal diffusion factor versus water molar fraction 
in water-ethanol mixture at 298 K and 1 atm. Simulation 
data: •, experimental data: □ ["HI . < [Ts|l . A [T^l . Dashed 
line is a guide to the eyes. 

In our simulations we have used a non-equilibrium al- 
gorithm pol l2l| , which permits to maintain a heat flow 
through the system with control on the average temper- 
ature and total momentum and energy conservation. As 
far as the molecular potentials are concerned, we have 
chosen simple models able to reproduce strong interac- 
tions such as hydrogen bonds. The following force fields 
were used for self-interactions: for water, the TIP4P 
model ; for methanol and ethanol, the OPLS model of 
Jorgensen |23|; for acetone, a model proposed by Wheeler 
and Rowley |24j: and, finally for DMSO, the force field 
derived by Luzar and Chandler [25j. Standard Lorentz- 
Bcrthclot mixing rules were used for cross interactions. 
The simulations were carried out on 800 molecules in 
water- alcohol systems and 500 in the other cases. The 
details of the simulations will be given elsewhere. 

The simulation results for the water-ethanol mixture 
at ambient conditions are compared with experimental 
data on fig. ^ It is important to notice that experi- 
mental data have been obtained by three different meth- 
ods 0, 0, 0| . The consistency between experimental 
results is remarkable and these data can be considered 
as reference data for this system. As can be seen, our 
simulation results are in very good agreement with ex- 
periments. This validates our methodological choices. A 
change of sign of ar is observed for x wa ter ~ 0.86. 

In fig. [3 we compare simulation data with experiments 
for the water-methanol system at 313 K and 1 atm. Ex- 
perimental data are those of Tichacek et al. [2(| at the 
same thermodynamic conditions. We observe a plateau 
at low water molar fraction followed by a change of sign of 
«t at £ wa ter ~ 0.8, in agreement with the experimen- 
tal data. Finally, predictions for the thermal diffusion 
factor versus water content in water-acetone and water- 
dimethyl-sulfoxide mixtures are shown in fig. [3 In both 
cases, a change of sign of ax is also observed. Therefore, 



FIG. 2: Thermal diffusion factor versus water molar fraction 
in water-methanol at 313 K and 1 atm: NEMD data (•) and 
experiments (O) from Tichacheck et al. [2t|. NEMD data for 
water-DMSO (▼) and water-acetone (A) at 298 K and 1 atm. 
Dashed lines are guide to the eyes. 

Despite the ability of MD simulations in reproducing 
the experimental behaviour, the question about what mi- 
croscopically causes the strong composition-dependence 
and, eventually, the change of sign, in ax remains still 
open. To shed some light on this question we have built a 
two-dimensional Ising model in a temperature gradient. 
In this model, the dynamic process consists of the swap 
of 2 particles between neighboring sites i and j, kept at 
different local temperatures T and Tj, according to an 
externally fixed temperature gradient. The presence of 
species A (or B) on a given lattice site i is represented 
by the state al = (1,0) (or (0, 1)) of the local spin. The 
site-site interactions arc limited to first nearest neigh- 
bors, and the lattice has no vacancies. The Hamiltonian 
is given by H = - >t <?j • J • (Tj, where J is a constant 
square matrix which defines the particle-particle inter- 
actions through its elements: saa, £bb, for the intra- 
species interactions, and eab = £ba for the inter-species. 
The necessary particle number conservation in the pro- 
cess implies a spin-exchange or Kawasaki dynamics |27|. 
The actual dynamic process is determined through the 
transition probability for one spin-exchange, which has 
been chosen to be 



W([a i ,a j ]^[a' i ,ar]) = 

min \ e -l( E *- Ei )/ kTi+ ( E 'i- E i)/ kT i] l 



(2) 



In this expression, Ei is the energy of the spin at the site 
i, that is 



(3) 



where nn stands for nearest neighbors. This transition 
probability satisfies the required detailed balance condi- 
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tion so that the appropriate thermodynamic equilibrium 
is reached if the temperature is uniform. 

In fig. |3| we present the results of simulations on the 
lattice model for different values of eab, with eaa = —2.0 
and ebb = —1.0 in dimensionless units, where T* = 2.55 
is the average reduced temperature. The figure clearly 
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FIG. 3: Thermal diffusion factor <xt for the mixture of par- 
ticles in function of the molar fraction of one component for 
several cross interaction parameters eab = —1.0 (•), —1.414 
(□), —2.0 (♦), —2.546 (A). Dashed lines are guide to the eyes. 

shows that the slope of the concentration dependence of 
o.t changes from positive to negative when eab varies 
from ebb to eaa- When \eab\ > \saa\, \sbb\, a change 
of sign of cxt with composition is in addition observed. 
This is again a strong suggestion that in dense systems 
the behaviour of olt with composition is severely influ- 
enced by interations between unlike species relative to 
interactions between like species. 

Obviously, the lattice model is too simplistic to draw 
conclusions, but it grasps an important aspect of the 
phenomenon. It can be therefore interesting to evalu- 
ate the role played by the unlike species interactions by 
doing direct molecular dynamics simulations on a sim- 
ple Lennard- Jones (LJ) mixture, which corresponds to a 
more realistic situation. We have thus computed cut us- 
ing the same MD methodology as previously described. 
Both species have the same mass and size (a = 3.405 A), 
but different energetic parameters (eaa = 0.996, Ebb = 
1.384 kJ mol -1 . Notice that the values for particle A 
correspond to argon). For the inter-species energetic pa- 
rameter, we have used a value of eab — 2.348 kJ mol -1 , 
much larger than parameters between like species. The 
thermal energy is RT = 0.998 kJ mol -1 . Simulations 
have been done in the dense liquid state and care have 
been taken that no demixion occured. A change of sign of 
cut is again observed. Effectively, for the following com- 
positions, xb = 0.1,0.5,0.9, the values of the thermal 
diffusion are, respectively, -1.97, -0.187, +1.679. The fact 
that systems of Lennard- Jones particles could exhibit a 
change of sign of cut with composition has remained un- 
noticed in previous works, partly because most of studies 
have made use of the standard Lorentz-Berthelot mixing 



rule, which place eab between eaa and Ebb- Obviously, 
the interesting physics underlying the behaviour of these 
simple systems under a temperature gradient deserves a 
deeper analysis and work is in progress. 

All these observations therefore suggest that in more 
realistic systems the chemical potentials at a given space 
point z are analogous to the energy Ei of a given particle 
at the i th site in the lattice model, in view of eq. ©. One 
may thus think the physical system as divided in small 
mesoscopic slices and the dynamics of matter transfer as 
being described by a Master equation. In analogy with 
the Ising model described above, one may argue that the 
transition probability for an exchange of 8ua moles of A 
at the slice located at z by jSub moles of B, initially at 
the slice z + 6 is proportional to 



W[z,z + 5] oc 



e -(p A (z)- lf , B (z))/kT(z) 
e -{^ A (z+S)- 1 ^ B (z+S))/kT(z+8) 



(4) 



The proportionality factor 7 is introduced due to the fact 
that the matter exchange is not necessarily equimolar. 
We have considered that the variation of volume in the 
switch of Sua moles by Sns should be negligible, so that 
vaStia + vbStib = and hence 7 = va/vb- In eq. (@J 
S should be considered as a characteristic microscopic 
length analogous to the mean free path in kinetic theory 
of gases. Respectively, ha and va together with [1b and 
Vb are the local chemical potentials and partial molar 
volume of species A and B. In steady state, driA/dt = 
and thus a probability balance indicates that 

{n A (z + 5)W[z + S,z}+ n A (z - 8)W[z - S, z] (5) 
-n A {z) (W[z, z-S} + W[z, z + 5})} = . 

Taking S as a length much smaller than the spatial vari- 
ation of the thermodynamic fields T(z) and ha(z), an 
expansion in powers of S reveals that, in a system in 
steady state with no net matter flow, the condition 



/ [i A - im 
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(6) 



should be satisfied. In the spirit of Onsager's linear irre- 
versible thermodynamics [28j, this result indicates that 
the effective thermodynamic force causing the whole ef- 
fect of species separation in a liquid in a temperature 
gradient is the left hand side of Thus, eq. @ states 
that in a system without matter flow, the variation of 
{ha — l^B)/kT with temperature is compensated by the 
variation of the same quantity with composition. This 
point of view permits to obtain a relationship between 
the temperature gradient and the composition gradient 
and, therefore, of the thermal diffusion factor, according 
to eq. Q, 
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vavb 



Ha/va - h B jvi 



xv a + (1 — x)vb x(dfiA/dx) 



(7) 
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where Ka and Kb are the partial molar enthalpies of the 
components A and B, respectively. In the derivation, 7 
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has been taken as constant. Surprisingly, this expres- 
sion is identical to that proposed by Kempers and very 
close to that of many authors that have used thermo- 
dynamic considerations (see [2^| and references therein). 
The accuracy of this expression has been tested on hy- 
drocarbon mixtures |scj . with a reasonable agreement. 
In associating fluids it is thus expected a change in the 
sign of the thermal diffusion factor with composition if 
Ha/va — h-B/vB changes sign at some particular point. 
This change of sign in our picture is caused by large inter- 
species interactions, larger than the intra-species inter- 
actions, which makes the dilute component to be more 
strongly bound than the concentrated component and 
then to have a more negative chemical potential. If the 
molar volumes are the same and, for instance, A is the di- 
lute component, any composition fluctuations will favor 
A to go downhill in temperature and accumulate prefer- 
ably on the cold side. This effect is reversed if we vary 
the concentration up to the point where B is the more 
dilute component: B is then more strongly bound than 
A and any composition fluctuations will favor an accu- 
mulation of B on the cold side of the box, thus reversing 
the sign of the Soret coefficient. If the molar volumes are 
different, any composition fluctuation involves different 
number of particles of both species present, to maintain 
the volume constant. Then it is essentially the value of 
the partial molar enthalpy per unit volume of the com- 
ponents, Ha,b /v a,b, which controls the sign of the sep- 
aration, as it has been already been pointed out in the 
literature |3l| . 

Equation © can be derived maximizing the appro- 
priate Legendre transform of the entropy, the so-called 



Massieu-Planck function S[1/T,P/T] H2, locally de- 
fined, and subject to the constraint of the constant vol- 
ume in the composition exchange. This is in fact the 
approach followed by Kempers. Under these conditions, 
the separation of the components is purely due to ther- 
modynamic effects, since (/j,a — 7A*b) I^bT is constant 
along the sample. However, the theory of irreversible 
processes also considers the existence of a contribution 
to the thermal diffusion due to purely kinetic origin [28| , 
which is completely ignored in the hypothesis underly- 
ing eqs. J2J) and (QJ . An explicit consideration of this 
kinetic effect would imply the right hand side of eq. © 
to be non-zero. Incorporating such a kinetic effects to 
the thermodynamic contributions related to the partial 
molar enthalpies is, to our believe, the main goal of a 
predictive theory for the Soret effect in describing both, 
chemical and kinetic effects, in dense systems as well as 
in dilute gases. 
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